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Abstract 

We are using a permanent-magnet MRI (Magnetic Resonance Imaging) with MRT spectrometer/imager to estimate 
the phase and void fraction distributions and flow patterns in gas-liquid two-phase flows. Air is being introduced at 
the bottom of the stagnant liquid column using accurate Harvard Apparatus programmable syringe pump PHD2000. 
Air flow rates were varied between 1 and 200 ml/min. The cylindrical non-conducting test tube of 2.5 cm diameter 
in which two-phase flow was measured was placed in a 10 cm diameter by 30 cm long solenoidal radio-frequency (rf) 
coil. 

The MRI is based on a 1900 lbs (865 kg) NdFeB low-field (0.267T or 2670 Gauss) permanent magnet protected by 
the Faraday cage. Very low fringe (stray) magnetic field exists which minimizes the rf disturbances and relaxes the 
safety considerations. The static magnetic field is very homogeneous due to the large (110 cm) return pole pieces 
that also minimize stray field. The operating proton frequency is about 11.3 MHz, which is outside the public FM 
rf range (it is in the HF range). Homogeneous region of 8 cm DSV exists with the homogeneity of 25 ppm p-p. 
Clear access of 13.6 cm exists in which the measuring sample, gradient coils and the rf probe-transmitter is placed. 
The resistivity of the pole pieces is larger than 1.2 /zOhm meter. The large magnet mass has a long temperature 
stabilization time and requires conditioned air with temperature control within 2-3 degrees Celsius for the stability 
of the static field. Intel Pentium PC computer with the I/O Texas Instrument TMS 320 C044 60 MHz DSP is used. 
F\ill functioning and flow imaging optimized Image-processing and spectroscopic analysis software is installed. Image 
processing software is based on IDL 5.6. 

The air is transparent to NMR rf signal and only precessing X H protons induce em/in the receiving coil. Accordingly, 
we measure gas phase distribution indirectly. The spatial signal is compared to a test tube filled with water to 
extract the gas-phase distribution. Fast temporal measurements of phase distribution in tube cross-section (s) and 
time-averaged void fraction for the 2D and 3D domains in case of a steady flow have been achieved. For void fraction 
imaging and data acquisition we used Hahn spin-echo technique with spin echo times varying between 6-10 ms. Due 
to the relatively fast echo sequencing and low bubble velocities good spatial and temporal resolution was achieved. 
The characteristic voxel size was less than 500 pm (0.5 mm). The spatial resolution can be improved down to 100 pm, 
in which case the temporal resolution would suffer somewhat. In Figure below, the measured vertically-integrated 
cross-section void fraction distribution in a vertical test tube with stagnant liquid at 30 ml/min air flow rate is shown 
(void fraction is on the vertical axis). 

MRI is a very promising technique for estimation of gas-liquid phase distributions, void fraction distributions, and 
flow patterns visualization in 2D or 3D domains for not too large flow velocities. A compromise has to be reached 
between the spatial and temporal resolution. MRI is especially well suited for 3D time-averaged volume- void-fraction 
measurements in case of steady two-phase flows. The beautiful thing with the MRI is that we can now add phase- 
contrast 3D velocity measurements by which independent closure of the mass and momentum balances can be 
obtained from the MRI measurements alone. 
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ABSTRACT 

We used a permanent-magnet MRI system to estimate the integral and spatially- and/or temporally-resolved void- 
fraction distributions and flow patterns in gas-liquid two-phase flows. Air was introduced at the bottom of the 
stagnant liquid column using an accurate and programmable syringe pump. Air flow rates were varied between 
1 and 200 ml/min. The cylindrical non-conducting test tube in which two-phase flow was measured was placed 
in a 2.67 kGauss MRI with MRT spectrometer/imager. Roughly linear relationship has been obtained for the 
integral void-fraction, obtained by volume-averaging of the spatially-resolved signals, and the air flow rate in 
upward direction. The time- aver aged spatially-resolved void fraction has also been obtained for the quasi-steady 
flow of air in a stagnant liquid column. No great accuracy is claimed as this was an exploratory proof-of-concept 
type of experiment. Preliminary results show that MRI a non-invasive and non-intrusive experimental technique 
can indeed provide a wealth of different qualitative and quantitative data and is especially well suited for averaged 
transport processes in adiabatic and diabatic multi-phase and/or multi-component flows. 
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INTRODUCTION 

The knowledge of the instantaneous spatially- resolved 
or time-averaged, with or without spatial information, 
phase and void fraction distributions is a crucial informa- 
tion in characterization of the multi-phase and/or multi- 
component systems. The flow patterns can be evaluated by 
the knowledge of the phase distributions and much about 
the flow dynamics in adiabatic or diabatic conditions could 
be learned. The degree of complexity in understanding 
mass, momentum, and energy transfer processes in multi- 
phase and/or multi- component flows in various geometries 
is enormous. Therefore, it is essential to provide as much 
accurate experimental insight and relevant data to better 
understand the complex multi- phase flow phenomena and 
phase interactions. 

Predictive CFD models rely on accurate and physically 
sound phenomenological relations, kinematic, and conser- 
vation principles for mass, momentum, and energy for each 
phase. CFD in gas-liquid flows have not yet reached the 
degree of sophistication required to compute phase distri- 
butions in complex geometries either in terrestrial or vari- 
able gravity conditions. The question is if there will ever be 
such a general two-phase CFD code that would account for 
all aspects of flow dynamics and energy transfer. Despite 
all the objective difficulties, there is a strong need to utilize 
two-phase gas- liquid flow systems in microgravity environ- 
ment due to their superior performance compared to the 
single-phase systems. The terrestrial two-phase flow sys- 
tems also need optimization that would be based on better 
understanding of the two-phase flow phenomena. 

Two phases of the two-phase gas-liquid system are liq- 


uid and gas or void phase. Void fraction is, by definition, 
the volume of void (gas) phase in a total volume occupied 
by both gas and liquid. Phase distribution is the instan- 
taneous spatially-resolved distribution of phases (gas or 
liquid) in a certain geometry. There is a need to measure 
void fraction on different scales for the whole volume do- 
main accurately, non-invasively, and non-intrusively. 

Current measuring techniques are usually intrusive 
and/or invasive, bulky, tedious, complex, or simply inad- 
equate to capture the dynamics of the flow. Some of the 
common non-invasive methods use a Photon-Attenuation 
technique (includes 7 -Ray and X-Ray radiation) that rely 
on the absorption law, I = Jo exp(— ^ a L), for the mono- 
chromatic collimated photon beam being attenuated over 
distance L. Obviously, such techniques can provide only 
line-averaged void fraction and in some instances area- 
averaged void-fraction by using multiple photon beams and 
detectors. The high-energy electron (/3-Ray) absorption, 
neutron scattering, and other particle radiation techniques 
are based on similar principles (( 2 , 13]). A large class of 
experimental techniques are placed in the category of opti- 
cal methods. Optical techniques require transparency and 
are plagued with multiple-scattering problems. Then there 
are probes based on the electrical properties of the multi- 
phase and/or multi-component mixtures, such as: conduc- 
tivity, capacitance, impedance, etc. Local probes such as 
micro-thermocouples, isokinetic probes, hot-wire and hot- 
film probes, etc. ([2, 13]), are also widely used but are ob- 
viously, both, intrusive and invasive, thereby affecting the 
measured flow. A major perturbation in the flow is caused 
by holdup technique (which rapidly seals off a section of 
a channel) which cannot be used for monitoring purposes 



but only for gross calibration of other techniques. The dis- 
cussion of shortcomings and advantages of each technique 
is beyond the scope of the present work. More information 
about particular techniques can be found in many refer- 
ences, such as [2, 13]. 

Magnetic Resonance Imaging (MRI) is a technique 
based on Nuclear Magnetic Resonance (NMR) and was in- 
troduced in the early 70’s. While the NMR is mostly used 
for spectroscopy and relaxometry, the MRI provides addi- 
tional dimensions and that is imaging or spatially-resolved 
signals to which many different imaging contrasts can be 
added. The inherent capabilities of an MRI are to measure 
the spatial distribution of the nuclear spins and as such 
needs a set of gradient coils that would, usually linearly, 
modulate the strong polarizing uniform field. The spatial 
resolution of the MRI can go down to 1-10 (xm while the 
temporal resolution can be typically on the order of 1-10 
milliseconds. Trade off between temporal and spatial reso- 
lution is required. The NMR is very insensitive technique 
with respect to SNR (Signal-to- Noise Ratio), but has a 
great advantage over all other tomographic techniques (X- 
Ray CT, etc.) because it provides unsurpassed level of 
contrasts. 

A review of MR imaging techniques as a tool to char- 
acterize various single- and two-phase flows is given in a 
number of excellent papers, such as, Majors, et al (1989) 
[16], Caprihan and Fukushima (1990) [4], Pope and Yao 
(1993) [17], Gladden (1994) [11], and Fukushima (1999) 
[10]. Most of the MRI investigations of two-phase flows 
is based on imaging of particulate flows and concentration 
measurements. None of the review papers reports any gas- 
liquid MRI measurements and/or research in the period 
from 1980 to 1999. The only two reported investigations 
using NMR (and not MRI) in gas- liquid flows will be dis- 
cussed now. 

The first paper known to us dealing with the use of 
NMR technique in gas-liquid void fraction measurements 
was done by Lynch and Segel [15] in 1977. The authors 
have used spectroscopic signal from the CW (Continuous 
Wave) NMR. Signal obtained is the Fourier Transform 
(FT) of the Free Induction Decay (FID), from the total 
volume that resulted in the volume- averaged integral- scale 
(sample-size) void fraction in the similar setup as ours. 
Air was introduced in a stagnant column of liquid. No 
spatial information on phase distributions or void fraction 
was provided as authors have used NMR instead of MRI. 
Nevertheless, it was an important contribution that proved 
feasibility in using NMR in characterization of gas-liquid 
flows. The experimental data show linear relationship for 
NMR signal intensity and void fraction versus air flow rate. 
The effect of the spin-lattice relaxation T\ on the signal de- 
cay was investigated as well. The authors have compared 
their integral void fraction measurements versus air flow 
rate with the homogeneous and Lockhart- Mar tinelli corre- 
lations. Authors did not discuss the problem of magnetic 
susceptibility. 

In historically the second paper on NMR use in gas- 
liquid flows, Abouelwafa and Kendall (1979) [1] have used 
CW (Continuous Wave) NMR spectroscopy to measure the 
signal intensity from the precessing spins ( L H) (protons). 


In this way they measured in situ the volume fractions 
and individual flow rates of two-component (water-oil or 
water-air) mixtures. The main goal of their research was 
to develop the flowmeter for the multi-phase and/or multi- 
component mixtures. At the time of our MRI experiments 
we were unaware of the papers by Lynch and Segel or 
Abouelwafa and Kendall. 

MRI has been successfully used in many investiga- 
tions of single- and two-phase fluid flow through porous 
media. We were motivated to extend the MRI tech- 
nique to gas-liquid flows in channels. In many ways the 
void fraction estimation is the inverse problem from the 
porosity /saturation measurements in porous media flows 

[5, 7, 8], 

THEORETICAL CONSIDERATIONS 

Two-Phase flow 

Two-Phase gas-liquid flow consists of liquid (holdup) 
and gaseous/ vapor (void) phase. Local or spatially- 
resolved phase distribution or fraction is designated as cq 
for liquid and a v for gas phase [13]. Since we are inter- 
ested primarily in the spatially- resolved or integral volume 
fraction of gas phase, we will use designation a = a v . In 
general, a — a(xi, t), i.e., local void fraction is a function of 
spatial and temporal coordinates. For the void phase, the 
smallest scale at which one can still assume continuum is 
on the order of 10 p,m. The liquid phase and the phase in- 
terface have smallest continuum length scales on the order 
of 10 nm , i.e., a factor of 1000 smaller than gaseous. 

We can now define a local instantaneous void fraction 
(or phase distribution function): 

. . [ 1 Xi £ void domain /1X 

a = a(x it t) = | o x< 3 vo id domain (1) 

If one could freeze the entire gas-liquid flow for an 
instant, we could perform volume-averaging or volume- 
smoothing process by increasing the size of the control vol- 
ume. That would result in a meso-scale volume 5V and the 
meso- scale void- fraction in a finite but small volume SV as 
a centroid around the point xf. 

6V 

e(zi,£) = J a(xi,t)dx? (2) 

o 

and is a function of time and spatial coordinates. For 
example, in the next instant different phase component 
can enter the Eulerian volume under consideration and/or 
phases are spatially distributed at each instant of time. 
The volume smoothing can go up to the integral scale or 
sample size resulting in an integral void-fraction: 

v 

<e(t)>= y r j e(xi,t)di cf (3) 

o 

which is the function of time only. This kind of informa- 
tion is accessible to NMR spectroscopy ([9, 10, 11, 15]). 
On the other hand, MRI can give additional information 
about the spatially-resolved void fraction, i.e., e(xi,t). In 
Fig. 1, we show schematically the thin cross-section slice 



NMR basics 
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Figure 1: The thin cross-section slice through the test tube 
with bubbly gas-liquid flow. 

through the tube with gas-liquid flow. Area-averaged void 
fraction can be easily estimated from such images. We 
can also define a meso- scale time r for the time- averaged 
void-fraction which will be the function of the spatial co- 
ordinate for the fixed Eulerian control volume: 

st 

e{xi,r) = \\m^ J a(xi,t)dt (4) 

o 

In the case of a quasi- steady gas- liquid flow: 

st 

e(xi) — lim — / a(xi,t)dt (5) 

St— *oo 01 J 
0 

which is the function of spatial coordinates only. Volume- 
averaging of this time-averaged void fraction results in 
both time and volume-averaged void fraction: 

v 

<£>= ^ J e(xi)6x f (6) 

0 

which one obtains in a steady gas- liquid flow by for ex- 
ample using holdup technique [13, 15]. All the smoothing 
(integration) techniques could be performed in terms of 
Lebesgue integral which exists even for the domains where 
the Riemann integral is not defined. In the case of steady 
gas-liquid flow < e >= < e >. 


The Nuclear Magnetic Resonance (NMR) phenomena is 
based on the non- vanishing quantum magnetic dipole mo- 
ment of certain nuclei [3, 9, 12, 14]. Proton or which is 
the member of the B ary on family of elementary particles 
with the mass of 939 MeV is also the most abundant nu- 
clei in nature). It also posses the strongest magnetic dipole 
moment which is fortunate since water is the most impor- 
tant technical and biological liquid fluid. The magnetic 
dipole moment tH — 1 Ju where Ji is the nuclear angular 
momentum or spin and according to the rules of quantum 
mechanics can only have discrete values. Nuclear magnetic 
spins are polarized longitudinally in an uniform magnetic 
field due to the dominant Zeeman energy-splitting interac- 
tion [3, 9, 12, 14]. For protons, which are also fermions , 
spin quantum number is ±1/2 and the orientation in a 
longitudinally polarizing uniform magnetic field Bo can be 
spin- up t (parallel) or spin-down j (anti-parallel). The par- 
allel orientation is the preferred lower energy level while the 
anti-parallel configuration represents higher energy level, 
with energy difference, A E = — hfo = —h'yBo — —huo. 
NMR and MRI experiments are notoriously insensitive at 
normal room temperatures because the thermal motion of 
the lattice will excite roughly half of the spins to ’’jump” 
to the higher energy level (anti-parallel |). The Maxwell- 
Boltzmann (MB) classical statistics can be used as a high- 
temperature asymptotic behavior of the quantum Fermi- 
Dirac (FD) or Bose-Einstein (BE) statistical distributions. 
The FD quantum distribution is correct distribution for 
fermions while the BE distribution is the correct quan- 
tum distribution for bosons. Fermions, such as proton, 
are elementary particles with half- integer quantum spin 
number for which Pauli’s exclusion principle applies and 
which have anti-symmetric wave function. Bosons are el- 
ementary particles with integer spin number (or zero) for 
which Pauli ’s exclusion principle does not apply and which 
have symmetric wave function. Photon, the carrier of the 
electromagnetic interactions, is a typical boson which pop- 
ulates the same quantum state with other photons, has 
zero rest mass, and zero quantum spin number. 


While the NMR can only detect volume- averaged (in- 
tegral) time- resolved void fraction < eft) >, the MRI can 
measure both, time- and spatially- resolved void fraction on 
the meso-scales for time and volume (also area- and line- 
averages). This inherent ability makes MRI more powerful 
than any other diagnostic’s technique. However, due to the 
Fourier imaging we can not achieve high temporal and high 
spatial resolution in the same time. Trade off is necessary 
for the particular case. However, fast imaging techniques 
such as EPI, FAST, snapshot-FLASH, DANTE, etc., [3] 
exist and will be explored. 

Additional advantage of the MRI is to measure the 
phase velocities it*’s using the phase-contrast techniques 
[4, 10, 11, 12, 16, 17]. The volumetric flow rate of each 
phase is: 

Qi = J UiOiidA i — l,v 

A 

Work is in progress to incorporate phase velocities and 
phase flowrates measurements with the spatially-resolved 
void fraction estimations which will be addressed in our 
subsequent publications. 


According to the MB classical distribution the popula- 
tion of parallel spins N a is related to the number of anti- 
parallel spins Np as: 

ft-*® (if) «. + »■>-" m 

Since kT » hujQ at room temperatures (T=300 K), the 
Taylor expansion of the exponential function leads to: 
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where B 0 is the uniform polarizing field in Tesla. Here, 
h — hj 2 tt is the angular Planck’s constant, wo = 7 B 0 
is the angular nuclear ( Larmor ) resonance frequency, 7 is 
the gyromagnetic ratio which is constant 42.58 • 2 n rad 
MHz Tesla -1 for 1 H, and k is the Boltzmann’s constant. 
At room temperature, the MB distribution show that ap- 
proximately only one in 100001 spins will actually result 
in net equilibrium magnetization, all other 10 5 canceling 
each other. Accordingly, in a population of 10 20 spins 
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only about 10 15 spins will result in a measurable equilib- 
rium magnetization which is why NMR/MRI is greatly 
insensitive- The relative population difference in high- 
temperature approximation results in: 


AN Na-Nfi _ h^Bo 
~W ~ N ~ 2 kT 


hf = fojjQ (8) 


and is a very small number at room temperatures in ther- 
mal equilibrium. This is because the lattice is much larger 
energy reservoir than quantum interactions at room tem- 
peratures. The NMR photons are in the radio-range, typ- 
ically / = 1 —>300 MHz, and thus of very low ionization 
energy. The 10 MHz radio- wave has wavelength of approx- 
imately 30 meters and its single photon energy is on the 
order of 10 -26 J or 10“ 7 eV. That explains why the MRI 
is so widespread in medical imaging. 


The quantum net magnetization in thermal equilibrium 
for *H is equivalent to the experimental Currie’s law: 

Mo = po B 0 oc ^ hu o « kT (9) 

4 k 1 1 

where po = N/6V is the spin density. Polarization in uni- 
form magnetic field will, however, take finite amount of 
time - called T\ or spin-lattice (or longitudinal) relaxation 
time. 


MRI basics 

The MRI was ’’born” in 1973 in seminal and indepen- 
dent papers by Mansfield and Lauterbur (recipients of the 
Nobel price in Medicine for 2003). MRI is based on NMR 
[3, 14] by adding ’’gradient coils” that can locally modulate 
the magnetic field and accordingly the local nuclear Lar- 
mor frequency to(xi) = 7 B(xi). This procedure is called 
frequency (or space) encoding [3, 9, 12] 

dB 

u(Xi) — CJq +7 ~Q^ Xi = + (10) 

where B 0 is a solenoidal field V • Bo = 0, and is also ir- 
rotational in the case there are no conductors inside the 
sample, i.e., V x Bo = 0. The difference of volume mag- 
netic susceptibility between two adjacent media (air /water 
interface) gives rise to distortions in the local magnetic 
field and thus to image artifacts. These distortions may 
be detrimental in imaging but can provide useful informa- 
tion about material inhomogeneity. The susceptibility x 
is a dimensionless tensor, which relates the contribution of 
the magnetic polarization M p of matter to the magnetic 
field H applied in vacuum, so that magnetic induction B 
in matter can be written as: 

B = ^ 0 (l + x)‘H = /z 0 (H + M p ) M p = yH 

We used Spin-Echo (SE) imaging [3, 12, 14] for the void 
fraction measurements. SE imaging in a gradient field re- 
sults in a spatially- resolved MRI intensity, i.e., the signal 
intensity from each voxel, which size can be changed, is 
distributed throughout the sample. The MRI spin-echo 
measurement yield the intensity of the NMR signal that 
is proportional to the spin density of the precessing pro- 
tons in each voxel. The amount of water in each voxel is 
directly proportional to the NMR signal intensity coming 
from that voxel. 
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Figure 2: ID x-axis cross-sectional spin-density distribu- 
tion in a cylindrical tube. Void fraction was averaged in 
a longitudinal z and y direction. Void fraction distrib- 
utions were obtained by subtracting two-phase flow sig- 
nal from the single-phase water signal (left). Difference in 
NMR/MRI signals for the single-phase and two- phase flow 
averaged void fraction distribution in a vertical tube with 
stagnant liquid at 30 ml/min air flow rate (right). 


The spatial resolution of MRI, 1/A Xj, is related to NMR 
absorption line width or spread in Larmor frequencies and 
is reasonably high in liquids: 


1 


7 

A Xi 
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The NMR signal intensity coming from the voxel filled 
only with water has maximum intensity, I T£ f- The 
NMR signal coming from the voxel on the meso-scale is 
Io(xi,t) < I re f . The local quasi- instantaneous void frac- 
tion is then: 


i(x u t) = Iref - 1 - < 1 (H) 

ire/ ire/ 

where the ratio Io(xi)/I re f represents the liquid hold-up. 
In Fig. 2, the cross-sectional longitudinally-averaged void 
fraction distribution is shown schematically. The esti- 
mation of the local void fraction in this way is straight- 
forward. However, there are a few problems. One is the 
inevitable attenuation of the signal and the inability to 
measure the spin-density at the time zero or at the mo- 
ment immediately following the tipping of the longitudinal 
magnetization. The first SE will yield smaller amplitude 
than the signal that would be coming from all in-phase pre- 
cessing nuclei at the time equal zero. The reason for this 
signal attenuation is in the irreversible, energy- conserving, 
spin- spin or T 2 (transverse) relaxation. The single SE in- 
tensity is therefore only a fraction of the original signal at 
time equal zero. Therefore, one has to repeat the SE exper- 
iment for different echo times and then fit and extrapolate 
the signal to time zero. This technique is called Multiple 
Spin-Echo imaging (MSE) [7, 8] It should not be confused 
with the classical CPMG ( Carr-Purcell-Meiboom-Gili ) or 
echo-train method [3, 9, 12] used to estimate T 2 relaxation. 

Sample was positioned in the rf (radio-frequency) coil 
inside the MRI and spins were polarized in a uniform Bo 
field after approximately STi’s. Then the equilibrium Ion- 



gitudinai magnetization is tipped by a weak Bi> rf it/ % 
pulse, at exactly the nuclear (Larmor) frequency, hence 
Magnetic Resonance, P recessing nuclei are inducing FID 
(Free-Induction-Decay) signal on the receiving coil. After 
time Ie/2 the ir-pulse in the transverse plane is applied to 
rephase the spins. The flip angle $ can be estimated from 
the tipping rf pulse Bi of duration r as: 9 = 7 B 1 T. 

Right after the ?r/2-pulse all components or the trans- 
verse magnetization are precessing with approximately the 
same phase. The net magnetization undergoes a coherent 
rotation. As a result of inevitable inhomogeneities in the 
local magnetic fields (T 2 relaxation), the coherence of rota- 
tion is lost due to the appearance of the different precession 
phases. Thus, the vector sum of all transverse magnetiza- 
tion components eventually vanishes. But as long as the 
precession frequency of each magnetization frequency re- 
mains unchanged, a simple permutation (time reversal) of 
the fast versus slow components will generate the reoccur- 
rence of the signal after another time t#/2. That permu- 
tation is achieved by the ?r-pulse. Accordingly, an echo of 
the FID signal will appear after the time corresponding to 
twice the pulse separation time [3]. This time is called 
Echo time in imaging. 

In the center of the Spin or Hahn echo, all magnetization 
components are refocused. Since the echo amplitude is 
smaller than the original amplitude at time equal zero it 
will look as if there are less 1 H nuclei in the voxel leading to 
void fraction over estimation. We are therefore using MSE 
[7, 8 ] to reconstruct the voxel’s proton density. 

The NMR SE intensity is proportional to the echo- time 
tE 1 

Ise = fa exp(— £ e/T 2 ) ( 12 ) 

where, fa is the unknown initial equilibrium magnetization 
(to be extrapolated) and obtained on the rf coil just after 
the ?r /2-pulse. Actually, there is a slight delay of 10-20 
us before the receiving signal is picked up in the rf coil 
due to the transmit- receive switch inertia. If, for example, 
we conduct two spin-echo experiments with variable echo- 
times Iei and t E 2 the corresponding spin-echo intensities 
will be Is ei and Ise2- Assuming a single-exponential de- 
cay, Eq. (12), one can estimate the transverse relaxation 
time 

T 2 = (Ie2 “ tEi) ( 13 ) 

*SE2 

and the initial magnetization (spin density) is then simply 
fa — Is bi expOjsi/Tb) Is ei < fa (14) 

For better accuracy and statistical confidence several SE 
tests with different echo-times were performed and the Iq 
was extrapolated by single-exponential fit using the Least- 
Squares fitting. 

The detected MRI signal in the k-space is the FT of the 
spatial spin density p(xi), also known as the 3D-imaging 
equation [ 12 ]: 

s(ki) = JJjcPxipix^e-^^ (15) 

Reconstructed image is obtained from the inverse 


Fourier Transform : 

p(xi) = j I J d 3 ki s(ki) e +i2 ’ rfel ' I < (16) 

Time-dependent component of the wave vector hi , which 
is also the Fourier- conjugate variable to the position vector 
Xu is: 

t 

ki{t ) = ~ [ Gi(r)dT (17) 

0 

In spin-waiy imaging used here the signal intensity de- 
pends on the repetition time tu and the echo time fr o f 
successive scans. The initial longitudinal magnetization 
M z ( r) depends on Ti relaxation time. Relaxation-time 
contrast can be applied to signal coming from each voxel 
location r, resulting in weighted signal: 

M (r. t Ri t E ) — M 0 (r) |l - exp - 

where Afo(r) is the thermodynamic equilibrium magneti- 
zation according to Currie law, Eq. (9). 

Three gradients are applied in the 2D FT spin-warp 
imaging [3, 4]. One is the slice-selection gradient G Zy then 
phase-encoding gradient G v for a fixed time interval, and fi- 
nally read-out gradient G x and acquisition of n data points. 
Phase- encoding gradient is then incremented and the whole 
procedure repeated m — 1 times after which a 2D FT of 
the n x m data matrix is performed for image computa- 
tion. To prevent signal aliasing by frequency folding during 
detection Nyquist sampling frequency must be used. 

EXPERIMENTAL SETUP 

Our MRI system, shown in Fig. 3, is based on a 1900 
lbs (865 kg) NdFeB low-field 0.267 Tesia (2670 Gauss) 
permanent-magnet protected by the Faraday cage. Very 
iow fringe (stray) magnetic field exists which minimizes 
the rf disturbances and relaxes the safety considerations. 
The static magnetic field is very homogeneous due to the 
four large (110 cm) return pole pieces that also minimize 
stray field. The health hazard limit of 5 Gauss is almost 
contained within the main body. The operating proton 
frequency is about 11.3 MHz, which is outside the public 
FM range (it is in the HF range). Homogeneous region 
of 8 cm DSV exists with the homogeneity of 25 ppm p~p. 
Clear access of 13.6 cm exists in which the measuring sam- 
ple, gradient coils and the rf probe- transmitter is placed. 
The resistivity of the pole pieces is larger than 1.2 fuOhm 
meter. The large magnet mass has a long temperature 
stabilization time and requires conditioned air with tem- 
perature control within 2 — 3 K for the stability of the sta- 
tic field, Intel Pentium PC computer with the I/O Texas 
Instrument TMS 320 C044 60 MHz DSP is used. Full func- 
tioning and flow imaging optimized Image-processing and 
spectroscopic analysis software is installed. Image process- 
ing software is IDL 5.6. Details about the NASA’s MRI 
test rig can be found in Daidzic ei al [7, 8 j. Detailed 
information on the NMR/MRI hardware can be found in 
[3, 9, 12). 

Ail of the experiments employed SE’s. The shortest echo 
time was 6 ms. For projection imaging, broadband (or 
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Figure 3: NASA GRC’s 0.267 Tesla MRI facility. 

hard) rf pulses of duration 120 ps were used. For slice se- 
lection, 5-lobe sine amplitude modulated pulses of duration 
1 ms were used. Two dimensional FT imaging was done 
with standard spin-warp sequences with imaging gradients 
of up to 1.9 G/cm. For quasi-steady flow's we repeated the 
measurements n « 2 — ► 4 times as the SNftoc yTL The 
NMR signal intensity is a j 3 BqPq/T [12]. In our MSE ex- 
periments, echo-time was changed from 6 ms up to several 
tens of milliseconds. Criterion for acceptable extrapolation 
was that I$ei had to be at least 80% of extrapolated Jo, 
and more than 6 different SE amplitudes must have been 
obtained. Minimum echo amplitudes used in fitting were 
larger than 20% of Iq. Water was doped with small concen- 
tration of CUSO4 para- magnetic salt to reduce longitudinal 
relaxation. 


Susceptibility mismatch across the air-water interface 
induces frequency shift. The image distortions around the 
interface will depend on the tensorial properties of the sus- 
ceptibility difference. For imaging with spin-warp a defi- 
nite distortion of the image around the interface will oc- 
cur. The frequency shifts caused by susceptibility mis- 
match can be resolved with spectroscopic imaging [3]. Air 
largely consisting of nuclei with zero magnetic dipole mo- 
ments is transparent to NMR signal and only processing 
1 H induce em/in the receiving coil by the law' of Faraday : 


emf 


-s/B' d * 

Coil Area 


Accordingly, we measure gas phase (void) distribution 
indirectly through the measurement of liquid holdup, Eq. 
(11). The spatial signal is compared to a test tube filled 
with water to extract the gas- phase distribution. Fast tem- 
poral measurements of phase distributions in cross-section 
of the tube and time- averaged void fraction for the 2D 
and 3D domains in case of a steady flow have been ob- 
tained. Due to the relatively fast echo sequencing and low 
bubble velocities good spatial and temporal resolution was 
achieved at lower flow rates. The characteristic voxel size 
was on the order of 500 pm (0.5 mm). The spatial resolu- 
tion can be improved down to 250 pm, in which case the 
temporal resolution would suffer. 


RESULTS AND DISCUSSION 

Nonconducting tube 20 cm long and 2.5 cm in diameter, 
shown in Fig. 5 (right), was placed in a 10 cm ID solenoidai 


rf coil of the MRI and an accurate syringe pump was de- 
livering air flow in upward direction at the rate varying 
from 1 to 200 ml/ min. In figure 4. the procedure for void 
fraction calculation from the spatial MRI signal is shown 
(see also Fig. 2). The two-phase flow signal is subtracted 
from water signal, voxel by voxel, resulting in void fraction 
for each meso- scale voxel, Eq. (refeg:mesoscaie). 




Figure 4: Difference in NMR/MRI signals for the single- 
phase and two-phase flow averaged void fraction distribu- 
tion in a vertical tube with stagnant liquid at 30 ml/ min 
air flow rate. 

The volume- and time-averaged void-fraction < c > 
can be estimated for not too large air flow rates by measur- 
ing the height of the stagnant liquid column Iq and height 
of the column when two-phase flow exists l\ > Iq, as shown 
in Fig. 5, resulting in: 

<e>=l-| AV = (h - 1 0 )^- (18) 

Integral void fraction calculated from the ratio of cob 
umn heights from Fig. 5, for air flow rate of 30 ml/min 
is 0.11 ± 0-01 by using Eq. (18). The void fraction ob- 
tained by averaging in two-different dimensions, Figs. 8 
and 9, show from the linear approximation at 30 mi/ min 
air flow rates void fractions very close to 0.11 which con- 
firms the visualization result. Unlike in [15], we do not 
need a (CCD) camera nor any optical access as the MR 
imaging capabilities are used for that purpose, although 
camera can be used for independent verification. 

Using the thin cylindrical slice perpendicular to 
tube’s longitudinal axis the cross-sectional holdup distri- 
bution is shown in Fig. 6 (left) for 30 ml/min air flow rate. 
Cross-sectional void fraction can be easily extracted from 
here. 


In Fig. 6 (right), a 3D image of gas-liquid Sow 
is shown. The blue-scale variation represents the void- 
fraction variation in different cross-sections. Since the cre- 
ation of such a 3D image with high spatial resolution takes 
long time (order of minute) it is necessary that the flow is 
quasi-steady. Each voxel’s void fraction is time-averaged 
according to Eq. (5). In Fig. 7, the projections of the 
3D image of two-phase flow void-fraction distributions are 




Figure 5: Side by side 2D images of the stagnant single- 
phase (water only) and the two-phase (water-air) flow at 
30 ml/min air flow rate. The gray scales are representa- 
tive of signal intensity (left). The schematic drawing of 
the test tube with stagnant liquid through which air was 
introduced is placed in the rf coil of MRI (right). 



Figure 8: Instantaneous 2D cross-section distribution of 
thin-slice longitudinally-averaged liquid holdup in a verti- 
cal tube with stagnant liquid at 30 ml/min air flow rate 
(left). 3D image for the particular slice of the spatially- 
resolved time-averaged void fraction for a vertical tube 
(right). 

shown. The flow occurs at high flow rate >150 ml/min. 
It is easy to see how the liquid film is contained at the 
wall forming annular flow conditions while the core of the 
cylindrical tube has high, and more or less uniform, void 
fraction distribution. It is expected that the gas-liquid flow 
in vertical flows with e > 0.20 will transition to slug and 
chum flow (bubbly annular transition). 

In Fig. 8, the variation of the void fraction with the 
flow rate (1 — 140 ml/min) is shown. The change is nearly 
linear for lower flow rates but for the higher low rates the 
MRI signal changes non-linearly and starts even decreasing 
for (> 180 ml/min). The reason for such behavior at high 
flow rates is the T\ relaxation effects and the short resident 
time of spins in voxels. The way to deal with this problem 
is to differentiate between stationary and moving spins. 
The integral void fraction was obtained by averaging void 
fraction along the tube’s length and then averaging along 
the y-axis resulting in the x-axis distribution. In the next 



Figure 7: 3D image for high air flow rate. 


figure, the distribution will be along y-axis, meaning that 
the results should be equivalent. However, discrepancies 
are possible due to misalignment, imaging sequence, gra- 
dient strength’s, and especially different concentrations of 
CuSC >4 para-magnetic salt. 
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Figure 8: Change of longitudinally (z-axis) and x-axis av- 
eraged void fraction as a function of air flow rate. Linear 
approximation shows reasonable trend. Large deviations 
appeared for higher flow rates. Errors estimated to ±15%. 


In Fig. 8, the variation of the void fraction with the 
air flow rate (1-200 ml/min) is shown. The change is nearly 
linear for lower flow rates but for the higher flow rates 
(> 140 ml/min) the MRI signal changes non-linearly. We 
already explained some of the reasons. Without the ve- 
locity information it is impossible to say just what should 
be the slope of the de/dm a gradient. The mass balance 
can not be closed by void fraction measurements alone. 
For very high flow rates the magnetic susceptibility ef- 
fects, aperture error, and Xi relaxation become significant. 
Better and faster imaging sequences will be used in the 
future to deal with such difficulties. Another problem in 
our imaging efforts is that the bubble formation was not 
steady but random. MRI was set to perform scanning in 
fixed intervals so that spatial aliasing could have occured. 
The bubble’s initial velocity and size were also variable and 
changing in upward direction. Bubbles experienced high- 
est acceleration at the formation site. This is also visible 
from Fig. 5 as the used MR imaging sequence did not 
characterize this region well 

CONCLUSIONS 

We demonstrated that MRI is a promising non-invasive 
and non-intrusive technique for experimental estimation 
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Figure 9: Change of longitudinally (z~axis) and y-axis av- 
eraged void fraction as a function of air flow rate. Linear 
approximation shows reasonable trend up to maximum air 
flow rates of 200 ml/min. Errors estimated to ±15%. 

of the spatially-resolved void-fraction. Also instantaneous 
or time-averaged void- fraction distributions in cross sec- 
tions and how patterns visualization in single- and multi- 
dimensional spatial domains is possible. We believe that 
void fractions up to 0.6 were estimated accurately with 
present imaging techniques. MR imaging has its inher- 
ent limitations in low sensitivity and is restricted to lower 
phase velocities. It needs non-conducting container walls 
and the image distortions due the magnetic susceptibility 
at the gas-liquid interface is present. A compromise can be 
reached between the spatial and temporal resolution. MRI 
is especially well suited for 3D/2D time-averaged spatially- 
resolved void fraction measurements in the case of steady 
two-phase flows or for the fast time-resolved ID void frac- 
tion distribution measurements. Gas-liquid flows can be 
imaged in adiabatic as well as in diabatic conditions en- 
abling the visualization of the void fraction change along 
the channel as the heat flux is applied. Another advantage 
of the MRI is that phase velocity measurements can be 
added by which independent closure of the mass and mo- 
mentum balances can be obtained from the MRI measure- 
ments alone. Additional improvements can be obtained by 
using different MRI contrast capabilities to visualize par- 
ticular transport phenomenon. The problem of magnetic 
susceptibility mismatch and image distortions in spin-warp 
imaging will be addressed by using spectroscopic imaging 
with phase encoding technique. More testing with differ- 
ent configurations, geometries, and flow parameters is re- 
quired. 
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